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Abstract - We continue the discussion on the interaction energy of the axially symmetric Hophons 
evaluated directly from the product anzsatz. The Hophons are given by the projection of Skyrme 
model solutions onto the coset space SU(2)/U(1). Our results show that if the separation between 
the constituents is not small, the product ansatz can be considered as a good approximation to 
the general pattern of Hophons interaction both in repulsive and attractive channel. 


Introduction. — Topological and non-topological 
' ^ solitons appear in many non-linear field models in various 

contexts. Since the appearance of the soliton solutions on 
pgthe theoretical scene in the 1970s, it has become evident 
that they play a prominent role in classical and quantum 
OO field theory. These spatially localized non-perturbative 
stable field configurations are natural in a wide variety of 
Cn physical systems [1] . 

The Faddeev-Skyrme model in d = 3 -I- 1 is a modi- 
^-nfied scalar 0(3)-sigma model with a quartic in derivatives 
^ term [2] . The structure of the Lagrangian of this model is 
• ^ similar with the original Skyrme model [3] , whose solitons 
are posited to model atomic nuclei, however the topologi- 
^ cal properties of the corresponding solitions, Hopfions and 
Skyrmions, are very different. It was shown that soliton 
solutions of the Faddeev-Skyrme model should be not just 
closed flux-tubes of the fields but knotted field configura¬ 
tions. 

The first explicit non-trivial Hopfion solutions were con¬ 
structed numerically by Battye and Sutcliffe m who 
found the trefoil knotted solution in the Faddeev-Skyrme 
model. Consequent analyses revealed a very rich struc¬ 
ture of the Hopfion states m- The subsequent develop¬ 
ment have revealed a plethora of such topological solutions 
with a non-trivial value of the Hopf invariant, which play 
a prominent role in the modern physics [5], chemistry [7] 
and biology [5]. A number of different models which de¬ 
scribe topologically stable knots associated with the first 


Hopf map S^ —>■ 5^ are known in different contexts. It 
was argued, for example, that a system of two coupled 
Bose condensates may support Hopfion-like solutions [9], 
or that glueball configurations in QCD may be treated as 
Hopfions. 

Note that most of the investigations of the Skyrmions 
and Hopfions mainly focus on the search for classical static 
solutions. Indeed, since in both models these configura¬ 
tions do not saturate the topological bound, the powerful 
technique of the moduli space approximation cannot be di¬ 
rectly applied to analyse the low-energy dynamics of the 
solitons. Therefore in order to investigate the process of 
interaction between these solutions one has to implement 
rather advanced numerical methods. 

Interestingly, the numerical simulations of the head-on 
collision of the charge one Skyrmions still reveal the well 
known pattern of the tt/2 scattering through the interme¬ 
diate axially-symmetric charge two Skyrmion m, which 
is typical for self-dual configurations like BPS monopoles 
[1]. However recent attempt to model the Hopfion dynam¬ 
ics m failed to find the channel of right-angle scattering 
in head-on collisions of the charge one solitons. 

Another approach to the problem of interaction between 
the stringlike solitons of the Faddeev-Skyrme model is to 
consider the asymptotic fields of the Hopfion of degree one, 
which corresponds to a doublet of orthogonal dipoles [H 
[13]. Investigating this limit Ward predicted existence of 
three attractive channels in the interaction of the Hopfions 
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with different orientation Ba¬ 
in his pioneering paper [3] Skyrme suggested to imple¬ 
ment so-called product ansatz to approximate a compos¬ 
ite configuration of well-separated individual Skyrmions. 
The ansatz is constructed by the multiplication of the 
Skyrmion matrix-valued fields. Note that besides the ra¬ 
tional map ansatz m it can be applied to produce an ini¬ 
tial multi-Skyrmion configuration for consequent numeri¬ 
cal calculations in a sector of given degree [16] . Evidently, 
the same approach can be used to model the configuration 
of well separated static Hopfions of degree one to approx¬ 
imate various multicomponent configurations. 

Recently we discussed the relation between the solu¬ 
tions of the Skyrme model of lower degree and the cor¬ 
responding axially symmetric Hopfions which is given by 
the projection onto the coset space SU{2)/U{1) [TT]. Us¬ 
ing this approach we made use of the product ansatz of 
two well-separated single Hopfions and confirmed that the 
product ansatz correctly reproduces the channels of inter¬ 
action between them. In this paper we briefly describe 
the relation between the solutions of the Skyrme model 
of lower degrees and the corresponding axially symmetric 
Hopfions which is given by the projection onto the coset 
space, adding here new and updated results of the numer¬ 
ical evaluation of the corresponding interaction energy of 
the Hopfions. 


rewritten in terms of coset space elements H, 

+ ^Tr{[a^H,a.H][a'^H,a'^H]}). 


(4) 


The difference between the Skyrmions and Hopfions is 
that in the latter case the dimensions of the domain space 
and the target space are not the same. The topological 
charge of the Hopfions, which meaning is the linking num¬ 
ber in the domain space |5], is not defined locally. 

There have been many investigations of the solutions 
of the model dTj) for higher degree Q |11[S1[T21[IH. Here 
we consider projections of general rational map Skyrmion 
ansatz, however numerical results of interaction potential 
is presented only for axially symmetric configurations of 
lower degrees Q = which are conventionally labeled as 
Ai.i and A 2.1 [5]. 

An rational map approximation to these solutions can 
be constructed via Hopf projection of the corresponding 
Skyrmion conhgurations HZ]. Recall that the rational map 
ansatz |18j is an approximation to the ground state solu¬ 
tion of the Skyrme model, which for baryon number B > \ 
takes the following form; 


Unir) = exp(2ih^foE(r)). (5) 


The model. — The Faddeev-Skyrme model in 3-I-I 
dimensions with metric (-I-, —, —, —) is defined by the La- 
grangian 

^ , (I) 

where (jf = {(p^, (fP, (p^) denotes a triplet of scalar real 
fields which satisfy the constraint |(^“P = 1. The finite 
energy conhgurations approach a constant value at spatial 
inhnity, which we choose to be (/)“(oo) = (0,0,1). For 
helds </>(x) with this property the domain of dehnition, i. e. 
3D Euclidian space, is equivalent to sphere and 0(x) 
dehnes the map ^ S^. It is well known that these maps 
are characterized by Hopf invariants Q = = Z, 

where the target space 5^ by construction is the coset 
space SU{2)/U{1). 

Any coset space element H can be projected from 
generic SU(2) group element U. In circular coordinate 
system the projection takes the form 

H = 2 = 2UToUt, (2) 


The unit vector nn is dehned in terms of a rational com¬ 
plex function R{z) = p(z)/q(z), where p{z) and q{z) are 
polynomials of complex variable 2 ; of degree at most IV, 
and p and q have no common roots. In Cartesian coordi¬ 
nates the components of n/j then can be written as 

Parametrizing the complex variable z = tan(0/2)e“'’ by 
polar and azimuthal angles 6 and ip we can Hnd explicit 
expression for any given rational function R. In particular, 
for baryon charge B = 1 we take R{z) = z, which gives 

fii = f = {(e”'"^ sin 9)/V2, cos 9, (—e“^ sin 9) / V2} . (7) 

We also implement the circular coordinates n±i = 
=F-^(h^ ± ih^) and ho = h^, where h® denotes the Carte¬ 
sian components of the unit vector (jS]). It is known that 
simple choice R{z) = z yields an exact solution of the 
model. 

For B = 2 Skyrmion the lowest energy rational map ap¬ 
proximation is given by the choice R{z) = or explicitly 
in circular coordinates 


where the Pauli matrices (ti,to,t_i) satisfy relation 


TaU = 7(-l)“(5a _f,l- 


V2 L 


I I I 

a b c 


(3) 


The symbol in square brackets denotes the Clebsch- 
Gordon coefficient. Then the Lagrangian (BD can be 


e ^“^sin^0 2cos0 — e^“^ sin^ 0 

•\/2(l -I- cos^ 0 ) ’ 1 + cos^ 9 ’ 1 / 2(1 + cos^ 0 ) 

For higher baryon numbers rational map approxima¬ 
tions are also well known [T], thought the correspon¬ 
dence between them and exact numerical results is get¬ 
ting worse as baryon number B increases. Note, that 
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in the topologically trivial sector B = 0 we may take 
no = {l/-\/ 2 , 0 , — 1 /-\/ 2 }, which we used to test the prod¬ 
uct ansatz approximation. 

Note that the Skyrme model can be consistently re¬ 
duced to Faddeev-Skyrme model by restricting SU(2) 
Lie algebra currents to coset representation ([2]). 

This means that Faddeev-Skyrme fields configurations 
with Hopf charge Q corresponds to Skyrme field with the 
baryon number B. Therefore the projection of rational 
map ansatz yields the rational map approximation of 
the Hopf charge of the same degree Q = N. Moreover, 
the usual profile function F{r) of Skyrme model, which is 
a monotonically decreasing function satisfying boundary 
conditions F{Q) = tt, F{oo) = 0, can be used. This pro¬ 
jection produces the Hopfion configuration of degree one 
with mass 1.232. 

As usual, we denote the Ai,i configuration as Hi which 
is a projection of the Skyrmion matrix valued field 

i.e. 

Hi(r) = 2Ui^..(r)roUjj^,(r), (9) 

where Ui{= 2 (r) is the usual spherically symmetric hedge¬ 
hog ansatz parametrised by rational map ([5]) with R{z) = 
z and hji is defined by eq. ([7]). It should be noted, that 
although the ansatz m for H = 1 is spherically symmet¬ 
ric, the corresponding Hopfion of degree Q = 1 does not 
possess the spherical symmetry. The projection breaks it 
down to axial symmetry m- 

The position curve of Hopfion is chosen to be the curve 
of the preimages of the point ( 0 , 0 ,— 1 ) which is the an¬ 
tipodal to the vacuum (0, 0,1). For the simplest Ai.i Hop¬ 
fion this is a circle of radius F{rc) = 7r/2, with numerical 
value Tc = 0.8763 in the x — y plane. Small deviations 
F{r) = F{rQ) + e then define the tube around the position 
curve where Ri 7 r /2 and (p = [0,27r). The same is true 
for R{z) = z^ or ^ 2,1 configuration, except that point on 
tube rotates twice when angle p changes from 0 to 2tt and 
the radius of circle being slightly large, = 1.299. 

For single Hopfion we can also rotate the points on the 
tube about the vertical X 3 axis by applying rotation trans¬ 
form via the SU{2) matrix 

H ^ D{a)liD{-a). (10) 

This global transformation corresponds to the symmetry 
of the Lagrangian (|1]). 

Let us now consider two Hopfions of arbitrary charges 
which are placed at the points R/2 and — R/2 and sep¬ 
arated by a distance i?, as shown in Figure „ There the 
polar angle 0 corresponds to the orientation of the Hop¬ 
fions with respect to the z-axis. Note that for i? = z the 
pattern of interaction between the charge one Hopfions 
is invariant with respect to the spacial rotations of the 
system around the z-axis by an azimuthal angle $. This 
additional symmetry does not appear for higher maps, e.g. 
for R = z^. 

First, we suppose that both separated Hopfions are of 
positive charges and they are in phase, i.e. for example. 



Fig. 1; Geometry of the system of two interacting Hopfions 
hi and ^ 2 , which are located at the points R/2 and —R/2 
(indicated by arrows), respectively. The dot position on hop¬ 
fion core illustrates the hopfion orientation (oposite phase case 
Aq = tt is shown here). 

rotation matrices D{a) in (fTHl) are identities for both Hop¬ 
fions. More strictly, the definition ”in phase” only means 
that the difference between the angles of rotations of in¬ 
dividual Hopfions is zero, Aa = 0. 

Then the system of two hopfions can be approximated 
by the product ansatz 

H^“=°(r) = 2Ufl.(r')Ufl..(r")roUl,4r")Uj,,(r'), (11) 

where r' = r -h R/2 and r" = r — R/2 and N denotes 
the sum of degrees of rational maps R' and R". Fields 
of both Hopfions in (EH) at the spacial boundary tend 
to the same asymptotics (0,0,1). Also note, that in the 
constituent system EH) of two Hopfions, contrary to the 
single Hopfion case, the transformation m of one of the 
Hopfions H do not leave the Lagrangian ((T|) invariant, it 
becomes a function of the relative phase Aa. 

In addition to the ansatz EH) we can as well consider 
the system of two separated Hophons, when one of them is 
relatively rotated by arbitrary angle a. Here, hovewer, we 
only restrict ourself to the relative phase Aa = tt, i.e., to 
the case when Hopfions have opposite phases. Taking one 
of rotation matrix D(0) to be identity matrix and the other 
D{tt) = 2to, we can express this system in terms of the 
matrices U/j' and U/j", thus the corresponding product 
ansatz is now different from EH: 

H^“="(r) = 8U^t.(r')roUfl,,(r")^oUlj4r")roU|j,(r'), 

( 12 ) 

The product ansatz approximations EH and EH en¬ 
sures the conservation of the total topological charge for 
any separation R and space orientation of the constituents. 
Note, however, that in the case of different rational maps 
R' ^ R", for example, z and z^, which are considered 
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below, their order in m and m is important. The 
different ordering yields different numerical results for 
small separation distances R, because the chiral angles 
for and TA//(r’") differ. This problem will be in¬ 

vestigated in more detail elsewhere. 

Substitution of product ansatzes m and m into La- 
grangian o yields energy densities of both configurations 
in terms of components of the position vectors r' and r” 
(cf Figure fl). 

Let us express these components via Hopfion’s position 
coordinates i?, 0, $ and the spherical coordinates r,9,ip, 
then the numerical integration of the corresponding lo¬ 
cal densities over the variables 'd and r yields the total 
energy (mass) of the system and its topological charge. 
The density of the latter quantity in circular coordinates 
is given by the trace formula 

Q{r',r'')=iV2i-ir+»\ ^ 

X Tr (V„ {Vr, (r')Ufl. (r"))U|j„ (r")U^, (r') 

X Vfc (Uh, (r')U^j. (r"))Ul,„ (r")Ulj, (F) 

We used the evaluation of the total topological charge 
of the configuration (for different values of the separation 
parameter R) as a correctness test of our numerical com¬ 
putations. The potential energy of interaction between 
two Hopfions of degree one and two is evaluated by sub¬ 
tracting of the corresponding masses of single Hopfions, 
i.e. rni^Ri^z) = 1.2314 and TO(_r"= 22 ) = 2.079503, from the 
integrated density (|3]) evaluated on the configurations (ITTl) 
and respectively. 



Fig. 2: The evaluated interaction energy of the Aq = 0 (in 
phase), product Q = 1 -|- 1 ansatz Hopfions as a function of the 
orientation parameters R and 0. 

Numerical results. — Evaluation of the total topo¬ 
logical charge and the energy of the product ansatz con¬ 
figuration requires numerical integration. In particularly, 
for each given set of fixed values of the orientation param¬ 
eters R, 0 and $, the integration of the energy density 


and Hopfion charge density (IT^ over three components of 
the Hopfion field yields, correspondingly, the strength of 
the interaction energy of the Hopfions and the topological 
charge of the configuration. 

We have performed calculations with different values of 
the parameters i?, 0 and $ for system of two product 
ansatz Q = 1 Hopfions [19] (figures [2] and |3|), system con¬ 
sisting oiQ = 1 and Q = 2 Hopfions (figures|l|-[7|), and for 
two Q = 2 Hopfions ffigures 151 - flUl). All figures demon¬ 
strate the integrated interaction energy as a function of 
the orientation parameters for the ”in phase” Hopfions 
and the Hopfions with oposite phases. Clearly, we can ex¬ 
pect our approximation of the interaction energy in the 
system of two separated Hopfions will be reliable only if 
the separation parameter R is larger than the sum of cores 
of the constituents Vc- 

Our results show that system of two Q = 1 Hopfions 
approximated via the product ansatzes m and (HU is in 
complete agreement with interaction pattern of the Hop¬ 
fions based on the simplified dipole-dipole approximation 
[13] . It should be noted that there was a minor mistake in 
our evaluation of the weight factors of the L 2 and L 4 terms 
in the effective Largangian presented previously in m- 
Fortunately, this mistake may affect the results only in 
the case of small values of the separation parameter i?, 
i.e. when the Hopfions cannot be considered as individual 
constituents, thus all predictions of m remain valid. 

In particulary, when the Hopfions are in phase and 
0 = 0, which corresponds to Channel A in [T5], there 
is a shallow (—0.05 at i? « 2.6) attractive window for 
separations R large than 1.8, as can be seen from Fig. [51 
Evidently, this attractive channel is very narrow because 
the potential of interaction quickly becomes repulsive as 
the value of 0 increases. 



Fig. 3: The interaction energy of the Aa = tv (opposite phases), 
product Q = 1 -I- 1 ansatz Hopfions as a function of the orien¬ 
tation parameters R and 0. 

When Hopfions are in side by side position, 0 = tt/2, 
the interaction potential is always repulsive as is shown in 
Fig. |5| The interaction strength for other orientations of 
the Hopfions is represented by a surface in Fig. |5| 
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Quite different pattern of interaction, however, occurs 
between the opposite phase Ai^i Hopfions, which is de¬ 
picted in Fig. [21 In contrast with Fig. [2] in the Channel A 
(0 = 0 ) the interaction is always repulsive for all values 
of the separation distance R. 

However, in the Channel B [T2] (0 = 7 r/ 2 ) the interac¬ 
tion energy has relatively large negative value (—0.266 at 
R= 1.6) at the separation of about two cores Tc = 0.8763 
(the interaction is attractive till i? « I.l) and then gradu¬ 
ally decreases to zero as separation between the Hopfions 
increases. Fig. [31 The repulsive behaviour changes to at¬ 
traction at 0 ~ tt/S, it approaches maximum as 0 = tt/2. 
For the Hopfions which are in opposite phases, this pat¬ 
tern is presented in Fig.jSjwhere we plotted the interaction 
energy as function of the orientation parameters R and 0. 
Qualitatively, the pattern of interaction between the Hop- 
Sons both in the Channel A and in the Channel B, is in a 
good agreement with results of full 3d numerical simula¬ 
tions of the Hopfions dynamics El- 

Unfortunately there is no similar 3d simulation data nei¬ 
ther for Ais, A 2 ,i configurations, Aip and A 2 ,i (figuresH] 
jT]), nor for two A 24 Hopfions ffiguresISl-ITUl) interactions. 
Below we consider the product ansatz approximation of 
these systems. 

It is clear, that in these cases the configuration are less 
symmetric and the interaction profiles are more involved, 
generally they depend on the value of the azimuthal an¬ 
gle 3>. In the case of Aip and ^ 2,1 Hopfion interaction, 
which in short will be denoted as 1 -I- 2, the shape of in¬ 
teraction energy isosurface for fixed values of the angles 
$ = 0 and $ = 7r/2 are shown in Figs [4] (in phase) and 
[5] (the case of opposite phases). We see that contrary to 
1 - 1-1 case strongest attraction (interaction energy mini¬ 
mum is —0.168 at R ~ 2.2) now is observed for Hopfions 
in phase (Fig.jT]), whereas for the oppositely oriented Hop- 
Sons (cf Fig [6]) there is only a shallow minimum (—0.014 
at i? ~ 2.4) for some narrow interval of values of the ori¬ 
entation angle 0 . 

In the case of interaction between two Hopfions ^ 2,1 (2-1- 
2 case), strongest attraction channel occurs again for the 
Hopfions with opposite phase, similar to the interaction 
pattern in the l-l-I case. The interaction energy minimum, 
however, becomes shallower with increasing of the Hopfion 
charges. In particular, in the 2-1-2 case the interaction 
energy approaches its minimal value —0.165 at i? ~ 1.4, 
compared to —0.266 at i? = 1.6 in the 1 - 1-1 case. 

In the case of the 1-1-2 system the interaction energy, 
in general, has nontrivial dependency on the value of the 
azimuthal angle $, as shown in Fig |S] In a particular 
case 0 = 0 (one Hopfion is above the other), the system 
possesses the axial symmetry and the dependency of the 
interaction potential on the value of $ is trivial. 

However for other values of polar angle0, for example 
for 0 = 7r/2 (the Hopfions are in the horizontal plane), 
this dependency can be evidently seen. Note that in the 
cases of the 1 - 1-1 and 2 - 1-2 systems there is no dependency 
of the interaction energy on the value of the azimuthal 


orientation angle $. 

In our consideration of the interaction between the Hop- 
Sons we mainly concentrate ourselves on the description of 
possible attractive channels. Note that all the possible re¬ 
pulsive channels are illustrated in Figs.BIlTUl thus we may 
also draw some conclusions about the repulsion strength 
between the Hopfions for different orientation angles and 
separation distance R. 

CD=0 




Fig. 4: The evaluated interaction energy of the Aq = 0 (in 
phase), Q = 1-1-2 product ansatz Hopfions is plotted as function 
of the orientation parameters R and 0 for fixed angles 4> = 0 
a) and 4> = 7r/2 b). 

*. — Conclusion 

We have investigated various interaction channels in the 
interaction of the axially symmetric Aij and ^ 2,1 Hop¬ 
fions. The product ansatz of Hopfions can be obtained via 
coset projection of the corresponding Skyrme field. This 
approximation preserves the topological charge in the en¬ 
tire interaction region. In particular, we analysed how the 
interaction energy depends on the orientation parameters, 
the separation R, the polar angle 0 and the azimuthal 
angle $. 

We have shown that this approach correctly reproduces 
both the repulsive and attractive interaction channels dis¬ 
cussed previously in the limit of the dipole-dipole interac¬ 
tions for Aip. Using the product ansatz we also were able 
to predict interaction pattern for pair of Ai,i and ^ 2 , 1 , 
and two ^ 2,1 Hopfions. In all cases the interaction has at¬ 
tractive channel for specific orientations and large enough 
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Fig. 5: The evaluated interaction energy of the Aq = 0 (in 
phase) ,(3 = 1 + 2 product ansatz Hopfions as a function of the 
orientation parameters R and "I? for hxed angles 0 = 0 a) and 
0 = 7r/2 b). 


separation distances R. 

Finally, let us note that the product ansatz can be ap¬ 
plied to construct a system of interacting Hopfions of even 
higher degrees and specihc spacial patterns. An approx¬ 
imation to the higher charge linked solitons, whose posi¬ 
tion curve consists of a few disjoint loops, like for example 
the configuration in the sector of degree <3 = 6, can be ob¬ 
tained as a multiple product of the projected matrices (0. 
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Fig. 7: The evaluated interaction energy of the Aa = tt (oppo¬ 
site phase), Q = 1 + 2 product ansatz Hopfions as a function 
of the orientation parameters R and il? for fixed angles 0 = 0 
a) and 0 = 7r/2 b). 


Fig. 9: The evaluated interaction energy of the Aa = 0 (in 
phase), Q = 2 + 2 product ansatz Hopfions as a function of the 
orientation parameters R and <1? for fixed angles 0 = 0 a) and 
0 = 7r/2 b). 


Oe[0,n/2] 



Fig. 8: The evaluated interaction energy of the Aa = 0 (in 
phase) a) and of the Aa = tt (opposite phase) b) of Q = 
2-1-2 product ansatz Hopfions as a function of the orientation 
parameters R and 0 for angles 4> € [0,7r/2]. 


0=0 



Fig. 10: The evaluated interaction energy of the Aa = tt (op¬ 
posite phase), Q = 2-1-2 product ansatz Hopfions as a function 
of the orientation parameters R and $ for fixed angles 0 = 0 
a) and 0 = 7r/2 b). 
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